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INTRODUCTION 

Spline functions are generally believed to have been introduced by I. J. 
Schoenberg in 1946. In the last decade, spline functions have attracted wide 
attention and the literature in this area has increased rapidly (approximately 350 
known journal articles, books or dissertations in the period 1961-1970). Unfor- 
tunately, many of the recent developments are not in a form which is convenient 
and accessible for application - oriented users. 

Spline functions, especially the cubic spline, have been a valuable addition 
to the fields of approximation theory and interpolation theory. Generally, they 
are much better than approximations which pass exactly through data points because 
they simultaneously approximate the function, its derivative and its integral. 

Actually, spline functions have been known for more than two hundred years 
and were first introduced in attempts to mathematically model the elastic curve. 
Shoenbeng did, however, introduce the name "spline functions" and was probably 
the first to systematically study the cubic spline and its generalizations and applicat- 
ions. There have been many generalizations of the cubic spline. But, for the most 
part, these generalizations had one common feature. They were generalizations of 
the cubic spline in terms of its properties, but not in terms of the original physical 
motivation. 

One mathematical basis for the development of spline functions is the Euler- 
Lagrange differential equation which arises by applying the techniques of variat- 
ional calculus to the problem of minimizing the integral 

b 

J = f y"7(i+y ,P f /p dx 

^a 

This integral represents the total strain energy in a relaxed, thin, elastic beam con- 
strained to pass, without buckling, through a prescribed set of points. Because of the 
complexety and non-linear nature of the fourth order Euler- Lag range differential 
equation which results from minimizing J, the exact problem has rarely been considered 
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Generally, there are a variety of ways to approach the study of spline 
functions. One approach is to grossly estimate the integrand in J and exactly 
solve the resulting problem. If the integrand in J is approximated by y" , the 
resulting problem lends itself to exact solution - the familiar cubic spline. 
Another approach is to investigate various approximations to the integrand in J 
and attempt to solve the resulting problems. This is the approach taken by the 
author. The results are described in the following pages. 

This report is a summary report on the progress of the author at the 
conclusion of approximately one half of the proposed project. Many questions 
remain unresolved. Most of the unresolved questions have received little or no 
attention - they were to have been considered in the latter stages of the project. 
The author does not feel that this report represents the report of a complete study 
of spline functions and has no intent of so implying. 
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MATHEMATICAL FOUNDATION 

Given a set of points ( x n , y^ ), n=o, 1,2, - — , N, where the x^, form 

a partition of the interval [a, b](i.e. a=x 0 <x<... x =b), we seek a family of 

1 n 

methods (functions) for fitting a curve exactly through these points. 

We desire further to have the functions from class C 5 [a,b], to have them 
satisfy a specified slope on curvature constrant at the end points, and to have 
them maintain the general global "shape" described by the specified points. 

We will restrict our investigation to a family of methods (functions) which 
are, in a sense to be described, approximations to the elastic curve. Assuming 
a desirable "shape" to be that assumed by a thin beam (spline) constrained to 
pass through the specified points and meet the end conditions, the resulting funct- 
ion is that which minimizes, relative to all admissable functions, the strain energy 
integral 



Our general method is to seek functions which minimize the integrals which result 
from replacing the exact intergrand in J by various approximations. The approx- 
imations which we will study are: 



2. y"V(l + 5y'/ / 2) 

3. y"*(l-5y'/$) 

4. y"7(l+5y’74) a 

5. y"7(l-5y'74? 

6. y" /(l + y‘ ) ^ (the exact integrand) 

For purposes of our study, we assume that y 1 is very small when compared to unity and 
certainly that y' < 2/5. The procedure will be to replace the intergrand in J by each 
approximation and apply methods from the calculus of variations to obtain a different- 
ial equation for the desired solution. 
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THE EULER-LAGRANGE EQUATIONS 


It is well known that a necessary condition for a function y to extremize an 
integral of the form I=J <f> (x / y / y' / y")dx is that y be a solution of the Euler- 
Lagrange differential equation 


3 / v 

— ^ / <>4>\ S4> 

dx S y9y/~ "3x ygy'J + ~%y 


= o 


V/hen viewed with the various integrands that we are studying^ it would appear that 
little success would be expected. For, in each case, the resulting differential 
equation is fourth order and non-linear. For example, if <t> = y" /(I + y 1 ) - the 
resulting differential equation is 


2y^(l + y' 2 f-5y' 1 [(y 1 '+4y*y ,,, )(l+y ,a )- 7 y ,r ’y ,,S ] = o 

i 

There is little improvement in simplicity by taking $ to be one of the other approx- 
imating integrands (excepts = y" z , which is rather trivial). 

Fortunately, in each of the cases we are considering, the approximating 
integrand does not explicitly contain terms involving y or x. This leads to a general 
method which permits two integrations of each of the fourth order equations. 

Since, in each case, there are no terms involving y, the Euler- Lag range 
equation can be written 



So, in any interval where y ,n is continuous, there is some constant 0 so that 
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ar(v) - V = “ 

where we have resorted to the usual subscript notation for partial derivatives. 
Now, 

IT = ‘V 7 ' + ^y* y " + V yl " 

But * y =° and * yl = jj_^ So, 

-^ +ayM =y " ^(v) + V r = i-( y "v) 


So, for some constant 0, 


4> + a y' + d = y" 4 > yll 

rb 

So, given the integral 1 = 1 <t> (y 1 , y") dx, where <t> has continuous first partial 

derivatives with respect to y' and y", the function y which extremizes I relative to 
all admissable functions (y'" continuous on {a, 6} except possibly at finitely many 
points a = x Q < x, < . . . <x^ = b and y € C" |a,bj ) must, in each interval 

[ X n-l, X i^ , n= 1,2, . . . N, satisfy the differential equation 


y" <t> ,, - <p = at y 1 + 0 
1 y" ^ n' n 

For the integrands which we are considering, this result can be simplified 
even further. In each of these cases, the integrand <t> has the property that 
y"</> y „ = 2<t> and the Euler-Lagrange equation becomes <t> = a^y' + 
Consequently, for the six integrands which we are considering, the respective 
differential equations are, in each interval jx^_^ x^j , n = 1,2,. . ., N: 

1. y" a = “ny 1 + ^n 
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2 . 

3. 

4. 

5. 

6 . 


y" 3 = (<v n y' + 0 n )(l+5y '%) 
y ,,f! = (at n y' + ^n)0-5y ,a /2) 
y" a = (o'ny'+ ^n)(l+5y'^) 2 
y"%(«n y'+ e n )/(\-5y fl /4f 
y ,,3 = (ar n y'+ /3n)(1+y ,2 fa 


For mathematical interest we will add two additional equations: 


2 (a) y 


» 2 = /3 n + Q - n y , + 5/3 n y'/2 


4 (a) y"°= /3n + « fl y' + 5 0nV' ? /2 + 5^/2 + 25 V*/ 16 


Note that the. right hand members of 4, 4(a), 2, 2(a) and 1 are polynomials in y' such 
that each is obtained from the previous by eliminating the highest degree term. 
Similarly, the denominators of the right hand members of 5, 3 and 1 are polynomials 
in y' such that each is obtained from the previous by eliminating the highest degree 
term. 
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SOLUTIONS 


Our attention is then focused to solving these eight non-linear, second 
order, differential equations. For convenience in notation during the discussion 
of solutions, we will assume that we are working in a specific interval |^ x n T x n ] 
and will delete the subscripts from the constants and Pn. 

The Equation y l|P = o-y* + [3 

By differentiation of each side of this equation, we obtain, where 
y" ^ 0, 2y'"=a.. Or, y^ 1V ^= 0. Consequently, every solution is a third degree 

ft P 

polynomial. Conversely, every third degree polynomial y = ax + bx + cx +d 

r> 

satisfies an equation of this type with a = 12a and p = 4b - I2ac. 

’ This equation is then an alternate approach to the well known cubic spline. 

Results relative to the cubic spline are extensive and well known and will not be 
reproduced here. . 

The Equation y” = (nr y' + ft)(1 + 5y'^2 ) 


For purposes of the general discussion, we assume that* / o and change 
variables by letting 


y = (2/50 ) (4tu - 5 Bx/6) 


The resulting equation is 


u)" = 4(a) ” e^) (to - e (a)'- e^ ) 


where 


j/ 

e^ = -5p/l2, e f =5 p/24 +i (5^/32) e^=e 
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This equation is the well known equation of the Weierstrass elliptic function. 
So, we have 


to' - t?(x+p) 

where the invariants of $ Qre 

= - 4 ( e^ e^ +e^ e^ + e^ e^ ) = 12(5p/24) P -4(5o ,! /32) 
g =4e e e =-8(5 p/24) f(5p/24) +5ry/32] 

9 1 r> 3 *- 

and A = g 3 - 27g <0 

r> r\ 

We now introduce the Weierstrossi “function ( £* (z ) = - )) # integrate the equation 

u)' = )P (x +P), and change to our original variables. From this we obtain 

y =A £ (x + B; g o , g^)+Cx + D 

where g and g are defined above, A = - 8/5a, and C = - p/3cr. . Direct substitution 

r> 3 

verifies that this is indeed a solution to the original equation. 

There is an interesting alternate approach to this particular equanon .in 
terms of Jacobi elliptic functions. Since y"= dy'/dx = y‘ dy'/dy, we can write the 
equation in the equivallent integral form 


x = s(2 /5) ly/ " | [(a u + p )(u% 2/5)] T/l, do 

and 

a y + P x = s(2/5) 1 /'’ J [(a u + p )/(u”+ 2/5)] ^ du 
where s = ± 1 = sgn(y"). By tetting u = a v/l a I, we obtain 

x= *j[(v-a)(v\b*]“ l/ *dv 
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and 


a y + Px = I rr I xj jjv-a)/(v 3 + b )] ^ dv 
where X = s(2/5)’/ 2 a, | a \^, a = - B / 1 a. I and b = (2/5)^ 

Now, we let A = (a +b ) 1 ^, k = £(A-a)/(2A)J ^ and 
v = A [(l-cn(u)/(1+cn(u)J + a, where cn(u) has parameter k. After some manipulation, we 

l/- j 


obtain 

and 


nr. 


= X A ' " 1 du 


y -t fix = I n I J |(I-cn(u)/(l+cn(u)J du 


Each of these integrals are known, and we have 


x + y * = X A" 1//s u 


and 


ay + 8x + 6* = I flt.1 X A ^ £u-?E(u) +2sn(u) dn(u)/(l +cn(u))J 
Here, Y*and 6* are constants of integration and E(u) is the fundamental elliptic 
integral of the second kind. We now let v =tu where t = s a/| a| =±1. Then, 
since t E(u) = E(v), t sn(u) = sn(v), dn(u) = dn(v), and cn(u) = cn(v), we have 

x + Y* =f?X5Ala 

and 

av + Px + 5 * = [2A I a \/5y [v-2E(v)+ 2 sn(v)dn(v)/( 1 + C n(v) )j 

With c ’= [2/(5A I a !)] and d 1 = (2A I a 1/5)^ ( so cc j = (5/2)1/=) 
we obtain 

cx + Y = v 


and 


Ctd y + ( P d-c) x + 6 = -2 [E(v) -sn(v) dn(v)/(l + cn (v))J 
This provides a functional relation between x and y. 
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The Equation y"~ = ,9 + a y'+5 Py'~/2 


. i/* 

First, we let X=(5f5/2 (which may be complex) and change dependent 
variables by letting y = a x + bM/ where a t~--ot/2 X and xb =( 2/5 -a This leads 

to the equation 

a " -x u,' => 

whose general solution is |i' = X sinh x (x + c )* Consequently, p, = cosh X (x + c) + c* 
and the desired solution is 

y =a x + b cosh x(x+c)+d 

In exactly the same manner, if we let h a = iy/2i and 

xb = (a 3 - 2/5) :/ we obtain 

y =a x + b cos X ( x + c ) + d 

as the general solution. These are of course equivallent solutions in the context of 
complex constants. 


The Equation y" = (a y* + ft )/(l -5y |g /2) 

We assume that a p / 0 and y'"< 2/5. Consequently, a y' + p *> 0. Now, 
we change variables by letting y' = a a m', where a = (5a /2) / ^ B , and obtain the 

equation 


0) 


,,r1 = ('"'-^/(a 3 -m' P ) 
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where b = -£/(a a ). From the restrictions y‘ < 2/5 and <xy' + P ■> O', we obtain the 
corresponding restrictions «>' -b”> 0 and a -<o 1 >0, Also, ab/0. 

Then, there are three cases that must be considered. They are: 

(i) b < -q < to ’< a, (ii) -a b <«>* <"a, and (iii) b = -a < (l) '< a. 

Case (i) b < -a<*» l< Q 

In this case we assume that y" (and hence u>") has no zeroes in the interval being 
considered and let s = + 1 be a factor to indicate the sign of m" in the interval . 

The differential equation is then 

sdm'= jj b -»« ')/(<•>' '-a*)J T/ ^dx 
or, ' • , / 

s ro ' d «> ' = £(b -«)')/(«>' - a )j' d 

So, we can write the two integrals as 

x=s £ |(a -u )( u +a)/(u - bjj ^'du 

and 

m= s J[(a -u)(u+a)/(u-bj ’ du t 


For simplicity, we write these as 


x = s ^[(a -u)(u+a)/( u - b)] ^ du 

“bx=s J[(a-u)(u + a)(u-fc)] ,//a du 


Now, we change variables of integration by employing Jacobi elliptic functions 
and letting 

sn(t) ~ £(a-u)/2a] ^ 
where the parameter k is given by k - [2a/(a - b)J ^ \ 


We then obtain 


= -2sk(2 a) 3 ^ J sn'(t) cn (t)dt 


n 


and 


« 


(u " bx = - 2 sk 1 (2cfl^ $ sn (t)dn (t)dt 

^ /P •**x 

Now, we simplify notation by taking X = -2 sk(2a ' andu = 
Then, we have 

Xx = Jsn (t)cn (t)dt 

and 

J4 ( (u -bx) =Jsn (t)cn (t)dn (t) c(t 


-2sk' , (2of /a . 


^ ^ ^ ^ *3 

Since dn (t) = 1 - k sn (t) and c n (t) - 1 - sn'XO, we obtain 

X x = A - A 

1 T> 4 

and 

u( m ~ bx) = A - (l + k')A +k^A 

^ 4 6 

where, again, A m = ^sn m (t) dt. 

Now, from well known reduction formulas we know that 
5k' A =4(l + k ? ) A - 3A + sn 3 (t) cn(t) dn(t) 

6 4 2 

Consequently 

P 4 

and 

5u(ni - bx) = 2A -(1 + k ) A f-sn S (t)cn(t)dn(t) 

' P 4 

Again, from known reduction formulas, 

3k A^ = 2(l + k 2 )A^- A Q + sn(t) cn(t) dn(t) 
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Then, we have 

3> k ? x = A 0 - (1 + 2k ) A - s n(t) cn(t) dn(t) 

and 

15uk ? («> - bx) = - (l + k*) A q -2(k 4 -k‘ + 1) A ^ 

-sn(t) cn(t) dnft) |7) k - 3k sn (t)J 

This result can be simplified further by writing it as 
3\ k ? x = A D - (1+ 2k ) A - sn(t) cn(t) dn(t) 

and 

15uk (m ~ bx)-3X k (1 + k ) x = -2(l + k ) A q -(k -1) A^+ 3k sn*(t) cn(t) dn(t) 

Now, k A^= A q - E(t) = t - E(t) where E(t) is the fundamental elliptic integral 

of the second kind. Thus, we have 

3\ k x = ( 1 + 2k ) E(t) - (1 + k ) t - k sn(t) cn(t) dn(t) 

and 

-15u k (•■> - bx) + 3>k (1+ k )x = (1 -k ) E(t) - (2k + 3k -1) t-3k sn (t) cn(t) dn(t) 
n.-i 

Now, '•> = (a a ) y, and constants of integration have not yet been added. So, we 
have the parametric solution 

Ay + Bx + C = (l-k") E(t) - (2k 4 + 3k 3 -l)t - 3k 4 sn®(t) cn(t) dn(t) 

and 

Dx + F = (1 + 2k 5 ) E(t) - (1 -i k S )t - k* sn(t) cn(t) dn (t) 
where a, p, C and F are integration constants, A = -15 uk (^ a ) , B = 3^ k (1 + k" ) 

+ 15ubk 4 , D^SXk 4 , k" = 2a/(a-b), r 1 =-2sk(2a) s/ ", 
u -, = -2sk “ T (2af^, b=-^aV), a =( 50 . 72 )"^ , and s = i 1 (= sgn? ( )"). 


13 


Perhaps the more useful form of this solution is 


and 


Case (ii) 


X x = Jsn (t) cn’ (t) dt 


(X + ub) x - u'h = k J sn 4 (t) cn P (t) dt 


-a < b < (0 1 < a 


As in the previous case, we arrive at the integrals 
x = s J [(a - u)(u +a)/(u - b)] ^ du 


and 


a) 


- bx = s^[(a-u)(u+a)(u-b)] 1 ^ du 


Now, we change variables by taking 
sn(t) =[(a-u)/(a-b)] 1 ^ 




where the parameter k is given by k = [(a-b)/2a] With this procedure we 
obtain 


and 


x = -2sk~(2a) ? j* sn"(t) dn“(t)dt 


bx = -2s k 4 (2a) J" sn‘ cn‘ dn dt 


or, with A ’= -2sk*'(2a) ,//u and u -1 — 2sk (2a)^= 2a k '^ \ 

>x = Jsn a (t)dn (t) dt 


and 


u(tt)-bx) =Jsn (t)dn‘"(t) cn“(t)dt 
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and 




u(nj-bx) = A -(l + k 3 )A +k n A 
z 4 o 

where 

A m =Jsn m (t)dt 

Now, from well known reduction formulas, we know that 

5k S A = 4(1 +k ) A 3A 0 +sn\t) cn(t)dn(t) 

6 4 z 

Consequently, 

X x = A k' A^ 

and 

5U(<u- bx) = 2A 2 -(1 + k^) A^ + sn 3 (t) cn(t)dn(t) 

Again, from known reduction formulas, 

3k 2 A 4 = 2(1 + k 2 ) A 2 - A q + sn(t) cn(t) dn(t) 

Then, we have 

3X x = Aq (1 " 2,k" ) A 2 “ sn(t) cn(t)dn(t) 

and 

15uk 2 (<u-bx) = (l+k")A 0 -2(l -k r> +k 4 )A 2 

- (1 +k“) sn(t) cn(t) dn(t) 

+ 3k%n 3 (t) cn(t) dn(t) 

This result can be simplified further by writing it as 

3Xx“ Aq+ (1 - 2k ) A^- sn(t) cn(t) dn(t) 

and 

1 5 1 - 1 k («» -bx) -3X (1 +k )x =-3(1 -k )A 2 t3k sn 3 (t)cn(t)dn(t) 
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p , 

Now, k A 2 = A q -E(t) where E(t) is the fundamental elliptic integral of the second 
kind. Thus, we have 

3Xk x = (l -k ? ) Aq- ( 1 - 2k )E(t)-k sn(t) cn(t)dn(t) 

and 

Xk (1 +k )x -5uk 4 (<•) - bx) = (1 -k ) A^+(l -k ) E(t) -k 4 sn^(t) cn(t) dn (t) 

This result can be inproved further by subtracting the first equation from the second 
to obtain 

3Xk x =(1 -k” ) Aq- (1 - 2k ) E(t) -k sn(t) cn(t)dn(t) 

and 

A k ? (k ? - 2) x - 5u k 4 (a) - bx) = (2 - 3k 3 ) E(t) +k P sn(t) cn(t) dn 3 (t) 

Now 

""T 

Ao = t, and we have not yet added constants of intergration. Also, =( a a ) 
So, we have the parametric solution 

Ay + Bx + C =(3k" - 2) E(t) -k sn(t) cn(t) dn^ (t) 

and 

Dx + F=(2k -l)E(t) + (l-k )t-k sn(t) cn(t) dn(t) 


4 :?\-T 

where <x, (?, C and F are constants to be determined, A = 5a k (o a ) 


B=-5Ubk 4 -X{ c a nc r ’-2), D = 3>k 3 , U P = 2a/(a-b), X” 1 
b = -p/(cr 2 a S ), a = (5 a " 2) ^ 8 ands= + l ( = sgn«>"). 


2 sk(2a) 


>/s 


4 /T* 

u = -2sk (2a) 


Perhaps the more useful form of this solution is 


Case (iii) 


b = -a < in ' ■ 


In this case, the differential equation becomes ">" n (a-fi>') = 1 , which is 
satisfied whenever 

frf =a - [+ 3 (x -n)/2] ^ 

That is, whenever 

(a a) 1 y = ax + & - 2^+ 3(k+y)/2l ^ A 


1 A 


where a , y and 6 are constants of integration, a=(5a /2) 7 and B - a a . 

Normally, these degenerate situations will not occur since there must be, 
in general, four undetermined constants in each solution interval. If the boundary 
conditions dictate the degenerate case, then the general case will in fact reduce 
to this case naturally. 


The Equation y" =( a y* <• p )(1 * 5y'/4) 

This equation is one of the more interesting of the equations we are 
studying. The general approach will be to obtain parametric representations 
for x and y and then eliminate the parameter to determine the functional relationship 
between x and y. 

We assume that a/ / o and note that rv y 1 + p > o with equal ity only at 
points where y" =o. First of all, we change variables by the substitution 

m = or y + p x. This leads to the equation 

fn" = > «)' («)' - 2p«)‘ + k 4 ) 2 
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where, for simplicity, we have let \ = 5/4 a and k = (f/ + 4cr. /5) 4 . We note 

o 4 

that the quadratic equation t-2P+k=0 has no real roots. Consequently, 
to ' ? Pm 1 + lc 4 has constant sign throughout the interval of consideration. We 

make the further assumption that ro" (and hence y") has no zeroes in the interval of 
consideration. 

s> dm'/dx = 2pu)‘ -i-k 4 ) 


A 


where s =+ 1 • Now, we may also write (since m" =du)'/dx = a> ' dm'/dtu ) 
s A (i , 1 dcn'/dti) = (m‘) ^ (d) 1 -2p(i)'-fk ) 


From these two forms, we conclude that 

du 


x = >s 


J 


^ ( n, - 2 p u + k ) 


and 


in - > S 


J 


di 


(n~ 2p u + |< ) 


where the parameter u. is actually <»'. Now, we change variables of integration by 

\/o 

letting v = (a) and obtain 


= 2*s j (v 4 - 2pv ' +k 4 ) dv 


and 


0) 


=; 2x s ^ v (y 4 - 2 p v +k) dv 


In this representation, the parameter v is actually (to/ ' = {<yy'+p) 1// > 
Fortunately,- each of these integrals can be integrated in closed form in terms of 
well known functions. After considerable effort, we arrive at the result that 
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X s 

— 

1 

/ V P+ 2k v cos (ft/2) + k ^ 

L 1 * -« 

/2kvsin(fj/2)\ 

a T C _ . a 

2k‘ 

2 co$(e/2) 

\ w ~ 2 k v cos (fl/2) + ky 

sin(ft/2) 



and 

to + d' = 


Xs 

TiT 


- 1 

/ v P + 2kv cos(n/2) + l< :, \ ] 

tan’/ 2l<v srn(9/2)\ 

2 cos (n/2) 

^ v' - 2 kvcos (P/^) + k ? y sin(G/2) 

l k ’- s jj 


z 1 / 

where cos = a fa / o < $ < IT , v = (u) 1 ) p , and c 1 and d 1 are constants of integration 
Consequent, we obtain 


0) 


+ d'+k°(x + c’ ) = sA sin (fi/2^tan 1 ^2kvsin(p/2)/(k -v 


and 


(to 


+ d‘ -k 3 ( x + c') = (> s/2k cos (o/2)) 


v T '~ 2kvcos(P*/2)+k 
v " + 2 k v cos ( 0/2) + k ^ 


Now, for breuity, we write 


and 


F (m, x) = exp [f ( (U , x )] 
G (d) , x) = tan [g (<*>, x)] 


where f( (U/X ) = [2kcos(fl/2)] [a, + d’ -k 9 (x +c‘)] / X s 

and 

g^/xH [k sin ( °/2)] [«, + d' +k"( x + c’ )] / x s 
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Then, we have 

v‘-2kvcos(e/5)+k r = ( v 3 +2kvcos(0/2)+k a ) F(«), x ) 

and 

2k v sin (0/2) = (k"-v') G ('», x) 

By simplifying each, we obtain 

v + 2k cos (fi/2) co th(f)v + k -0 (*) 

and 

v + 2k sin ( 0/2) cot(g) v - k^ = 0 (**) 

If we add (*) and (**), we obtain 

v + v k [cos (0/2) coth (f) + sin (0/2) cot (g)] = 0 
so (since v=o is not a solution) 

v = -k [cos (0/2) coth (f) + si n (0/2) cot (g)] 

If we subtract (*) from (**) we obtain 

v = -k/ jcos(0/2) coth (f) - sin (0/2) cot (g)J 

Since vA >0 an< ^ cos ( p /2) > o, we conclude that coth (f) < o (and hence f < o ) 
Equating the latter two expressions, we eliminate v and see that 

{cos (0/2) coth (f)j - {sin (0/2) cot (g)] =1 

Consequently (since coth (f) < o and cos (q/2)> o ) a necessary condition is 
- cos (0/2) coth (f) = [l+ sin (0/2) cot (g 2 ] l A 
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where 


f = 2k cos (0/2) [ay + (j 3 -k 3 )x+c] /ys, 

g=ksin(fl/2) [ay + (p+k*)x + d] /\s, 

sin(0/2)+[(k r '-P)/2kT /:? , 

cos ( 9 / 2 ) =[{k"+p)/2k 

1 / 

k = (p 1 - 4a°/5 ) 4 
X = 4 n -/5 , and s = sgn (y") = f 1 

The constants nr. , P, c , and d are arbitrary constants to be determined by boundry 
conditions. 

We note that the above result gives a necessary condition. For 
sufficiency, assume that x and y are related as above. Define a function p. (x,y ) 
by 

a ( x, y ) = sinh 1 [sin(0/2)/t an g] 

Then, sinh |i=sin (0/2)/tan g and cosh u =- cos ( 0/2 )/tanW f. 

_ u 

Now, we let v = k e >0 and determine that 
tang = 2k v sin (0/2 )/(k -v ) 

and 

tanh f = - 2kvcos( R /2)/(k r? +v°) 

Now, we solve for f and g and differentiate with respect to v. The result is 
' 7 y + (p + k*)x=2Xs(v-tk )/(v-2k cos 0 v + k 4 ) 

and 

rr y + (p-k ) x = 2x s( v -k ')/( v 4 - 2k cos 0 v +k ) 

Then, 

nry + px = 2Xsv r ’/(v^- 2 R v'+k 4 ) 
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and 

x = 2Xs/(v 4 - 2 fi v 3 +k 4 ) 

Then, we have 

*y'+M v P 


and 

y" = (dy 1 / dv)/x = v(v - 2p v +k )/a>. s 


So, 

(«y' + p + 5y ,t> /4)~ 5 v(v 4 -2r V*+ k 4 )/4a P 
Since a \ = 4 <* a /5, 

sy"=5v(v -2p v S+ k 4 )/4r, 2 

Consequently, 

y"=s(«y' + P) l/a (l+5y' a /4) 


The Equation y""= ft + ?. y' + 5Py ,r yfe + 5g.y ,fl /2+ 25By 4 /l6 

As in other cases, we assume a. ft ^0 and let s = sgn(y") = + 1 . 
y" = dy'/dx = y 1 dy'/dy, we obtain the integrals 

x = s f(a _ u* + 4a ' u" + 6 a u + 4a u + a ) 1 ^ du 
J' ° 1 r* ? 4 ' 


and 


= s .f u ( a o u4 


+ 4a 


* 

u + 6a 


A , 

+ 4a u + a ) du 

.*» 4. 


Since 
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where a c = 25 8/16, a = 5 a. /8, a =5 8/12, a =n/4anda = ft 

The parameter u is, in fact, y*. Now, we assume that y" = 0 when y' = b D 
(i.e. b Q is a root of the quritic); and let u = b Q + u>” 1 in the integrals. After 
considerable algebraic manipulation, we obtain 



and 


where 


x = -s $ (4t ? - g t-g ) % dt 

* o 

(y-b 0 x)/A i = -s $ (t-A /2) 1 (4t -g t -g^) Y° dt 
9n = 3Ag - 4A A = a 0 a -4a a + 3a 

~ a 1 .7 ° 4 17 7 
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and 


If *7 *'! •-» 

g ~ 2A A A - A - A A = a n a a +2a a a -a - aa -a a 

? 1 2 .7 T 4 u a 4 1 S :■ O o l 4 

Note that g and g are, in fact, the fundamental invariants of the original quartic 




f(y'). In terms of n . and 8, 


9„ =(25BV>2) -(5a7 8) 


9, = (58/6)" - (5 8/6)(5 a n /8)(1 1/1 6) 


Once more we change variables by introducing the Weierstrass Pe function 
$ (z; 9„,9 ) = )?(z) and letting t = $(z). This leads to the integrals 


x = s$dz 


(y- tv/ )/A i = s $ ( $( z) - A a/2)"^ 


The second integral is the Weierstrass elliptic integral of the second kind. With b 
chosen so that A = 2 i?(b), we have 


x + y = sz 


R + (y-b o x)/A j = [y (b)] 1 |xn [f7 (z-b)/ T(z+b)] + 2z r (b) J 
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Here rr and C are, respectively, the sigma and zeta functions of Weierstrass. Because 
the sigma function is an odd function (i.e. <r (u) = - 17 (u) ) we obtain, for either 
value of s, 

+ (y"^ x )/ A ! = L^'( b J ^ [ a (x + V-b)/rr(x + y + b)J + 2 C(b)(x+v)| 

and y has been expressed as a function of x and four constants (a/ Y, and 5) 
integration. The sigma and zeta functions are formed using the invariants g^ and g^„ 

The appearance of the numbers b and b Q in the solution will present great difficulty 
in fitting boundary conditions as will the dependence of g ? and g^ on the boundary 
conditions. 

For sufficiency, we note that 


(y'-b 0 ) =={*/# '(b)] j[rr '(x + Y-b)/a(x l Y-b)] - [a'(x + y +b)/a(x + y + b)] + 2a «(b)/cr(b) 
= A/ [5P(X+ V)-Sf*(b) j =A /[ S (X+V)-A/^ 


So, 


y M = -A x ft '(x + Y)/ [g>(x + y)“ A o / 2] 


Consequently, we have 

y "*=( y’-kj [»'(k + y)]7a’ 

But, 

(y (x + y)] “ = 4 y (x + Y) - g ( x + y) - g^ 
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So, 


y " ~ \ (y ' ' b o) 4 [ 4 ^{x + v) - # ( x+ y) _ g J 

However, 

j?(x + Y) = (A c /2)+A /(y'-bo) 

and we substitute to obtain 

■ / 

/“"-A/WI 4[A/(y'-b 0 )]’ + 6A r , [a /(y'-bof } 

* PA’- S^fA/fy'-bo)] + 4(A/2f-(g ? A/2 ) - ^ 

Now, 3 A -g =4 A A and 4(A /2) - (g A /2) -g =A A. 

2 3 1 •> V 2 2 fl 1 4 

So, we have 

y" r = A^ (y'-bof + 4A^(y , -b 0 f + 6AJy , -b 0 f + 4A i (y*-b 0 ) 

Expansion and simplification gives the original differential equation (since f (b D ) = 0). 


The Equation y” = (a v' + fl)/( 1 -5 y’ ^4) 

This equation lends itself well to a parametric solution. However, no 
method has been found to eliminate the parameter and finda functional relation 
between y and x. 

As in previous cases, we assume nr. p f o and that y" has no zeroes in 
the interval of consideration. 
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Then, 

y« = sjrr y'+ft) 1 ^! -5y' /A) 

where s = -t 1 = sgn ( y" ) 

We now change variables by letting y + ^x. Then, 

U) " = 4 a , s(u) , ) T / s /(4a -5f^ + 10f(ju' - 5a>' ) 

Then, we can write 

x =- (4 a ' , s)" l |[5,r + 10jp p -(4^ -5f )] iTA|u 
and, since du)‘/dx ~ (d tu '/cl(u) u), /i 

(o=_(4a q s)' 1 J [5^ ~ 10 R ii -(4 a - 5 fiP )j u 1 / 2 du 


where y, i s actually to’ 

Now, we change variables in the integrals by taking t = and 

integrate, to obtain the parametric solution. 


and 


ai 


where the parameter t is actually (ay' + ji) 


+ V = -(2^s) 1 [V’-10 R t y /3 - ^a P )tj 

+ a = - ( 2 a? s)~ 1 [5t 7 A-2Pt B -(4a P -5p P )t ? /i] 

1/2 


We note that, regardless of any knowledge of the significance of 
the parameter t , t x = <u . From this relationship we can easily verify that these 
parametric solutions are indeed sufficient to satisfy the original differential 
equation (assuming that t is positive). The significance of this parameter will be 
of considerable interest when we consider the boundary conditions and the 
problems of joining solutions over adjacent intervals. 
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The Equation y" =(ay' + 8 )(1 + y 1 ) 


V 3 


In this case we assume that, in the interval under consideration, 
rtfif o . We effect a rotation of coordinates by making the change of variables 


(i = ( a x-Ry)/(n”+fr ! ) ? 
v= (f x +a. y)/(n^ + p ) 


Then, l* =du/dv= (a. - py')/( a y* +p). Interestingly, this gives 
y' =( a ~p M,)/( a >* + f). So, 

y" = \r (a + P)V( a U + P) , 

nr y' + f = ( n ''+f)/( r <V- + p) " ° / 

and 

1 + y'" = (n- + p n ) ( I +{? + pf 

Hence, the differential equation becomes 


fl” = 

Now, we let 

u = 2 d)/(a. + P ) 


v = 2T/(« a +P ? ) l/4 

Then, A - d(o/dT = d|i/d v andu = d°u/d y = (d^u/cTT ) ('Af and we have 

the differential equation 


of =4(1 +m ? ) 5 / s 
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By direct substitution, we discover that the equation is satisfied if 

m = (m + y) “ 1 

where Y' s a constant of integration. 

Now, we again change variables by introducing the Jacobi elliptic 
functions and letting 

(0) + v) = cn (s,k); k = ^ 

We then obtain the differential equation 

(dT/ds) = k cn 4 (s,k) 

So, 

T = + k j" cn (s, k)ds = * k j" ( 2d n (s, k ) - 1 )d s 
Hence, T+Y=+k(2E(s, k ) - s) 

where E (s,k) is the fundamental elliptic integral of the second kind. 

By writing a = n/2(a and b = f/2(a.% we can write this solution in 

the form 

(bx+ay+v) = k (2E(s)-s) 
s 

where E(s) = j" dn^(t,k) dt, cr<‘ (s,k) =(ax - by + a ) ? and k = 2 

o , z 

To determine sufficiency, we note that d (2E(s) -s)/ds = cn (s,k) so that 

(b + a y') (bx +a y + s) =k 3 (2 E (s) - s) cn ? (s,k) (ds/dx) 

But, -cn (s,k) sn (s, k) dn (s,k) (ds/dx) = (ax - by + y ) (a - b y'). 
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So, 

(b + ay'f =k a cn* (s,x) (ax -by + Y) ? ( a -by 1 ) /cn (s,k)sn (s,k)dn (s,k) 

= k ? cn 4 (s,x) (a -by') /(I - cn^ (s, k) ) (1 -k sn (s,k) 
c 4 n(s,x)(a-by')V(1 -cn 4 (s,k)) 


Or, 


(a -by‘)V(b + ay') =L 4 - 1 

where L=ax-by+y • From this equation, we obtain algebraically 

(a%b r )(l+y ,7, )=L* 4 (b + ay') P 

and by differentiation 

(a +b )y" =2L~ S (b + ay')* 

Consequently, 

y" (1 + y‘ ' = 4 (a +b )^ (a y* + b) = ay* + B 

Several remarkable results can be observed by rewriting the above solution as: 

ax -by + y = f( t) = t cn(t,k) 

bx+ay + Y - kg ( t ) = k I cn a (u/k)du 
o^ 

where, g=f 2 (0, k=2 ^ , a = rr./2 (a + p S ) ^ and b= ^/2(a + P ) ^ . 

We again check to verify solution of the differential equation by observing that 
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ax -by = f 
bx +a y = kf' 

r=k* o-f 4 ) 

? + f ? =o 

x = (af + bkr)/(a"' + b*) 
y = (-bf + akf )/(a 2 +b S ) 

ay'+fl =kf (rr a + pb)/(bkf -t-a f) 

bk f + a f> 0 

y , *-kf(a%b' , ) 2 (bkf%af)”^ 

S gn(f)=sgn(y") 

1 +y ,r = k P (a a +b S )(bkf ? +af)“'’ 


We note especially that the curvature K is given by n = y"/(l + y ,rV ) Va = 


i/s 


.'fa 


2(a'+b < ) ' f = 2(a +b ) ' ” (a x - b y + v). Of more interest, and of considerable 
importance, is the significance of the parameter t. IF s is are length, then 
s =x +y = (f 5 + k‘ f 4 )/(a ?? + b ) = k /(a + b ). That is, the parameter t is, in 
each interval, a constant multiple of arc length. We have then derived a 
parametric representation of the true elastic curve with arc length as the parameter. 
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MINIMIZATION OF INTEGRALS OF POWERS OF CURVATURE 


The solution process for the previous equation, which arises from minimizing^ |<’ds 
(where k is curvature as a function of s), suggests a method of approaching a much more 
general problem. Though not directly a part of the present study, this generalization 
will be presented. 

Suppose that we wish to find a function y £C [a,b] which extremize the integral 
Jk n ds where k is curvature, n is a positive integer and n 5^ 1 . From our earlier work, 

we see that the Euler-Lagrange equation is 

n -i (3n-l)/2 

(y") = (n - 1 ) (a-y'+ 0) ( I+y 1 *’) 

Now, we change variables by making the substitution 


With u = du/dv, we obtain 


a/j = o x - Py 

a v = p x +£* y 

a = (n-I) ^ (a ? + p^) 


n n „ ( 3 n-D/2 

(o) * (-D (i +a ) 


(n-I)/2n 


By direct substitution, we find that this equation is satisfied whenever 

2n/(l -n) 

u'= |(n-l)p/n + Y*J -1 


-i / 2 

Now, we put k = 2 and 


Then, we obtain 


|(n-l)u/n + Y*J 3 = jcn 4 (t,k)j 


(n-l)/n 


(dv/dtf = 4k C [cn 2 (t,k)' 


(3n -2)/n 


For sufficiency, we write 


"x-py+Y - Ar/(n-l)[f(t)] 


(n-l)/n 


t (3n-2)/2n 

p x+oty + = 2kBj[f(t)] dt 
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where f(t) = cn (t / k) / k = 2 ? , A° = (n-l)(° (P + 8 3 ) and B = i A 

We observe that 

"x- By -Aff" 1//n 

(3n-2)/2n 

fx + 'Tfy = 2kBf 

f* =4k n f(l-f ) 
f = 2k (1 -3 f ) 

(a% f) = a. Aft 2kpBf^ 

(« r ' + f = -pAf i-2k 3 B f' //s 

ay'+p = 2kB(n r + p^f’^A^kpBf^) 

1+ y *** =4(n' , + P 5 )B n k n f/(" Af + 2kPBf 

(3n + 2)/2n . 

y" = (2kf AB(n z +ff f /(a A f+2k0B f y ’f 

,(3n-l)/2 

Direct substitution shows that indeed, (n-l)(y”) =(ar y' +p)( 1 + y‘ ) 
Alternately, we can write * 


« X - P y + Y = [A n /(n- 1 )] [f( t) J ^ n 


and 

l , _ (n-i.)/h n . 

px + ay+*=A| [f(u)] <l u 

o . 

where F(t) = cos^t and A =(n-l)( a +0 ) ^ ^ 

We observe that 


n.i- f!y = l/2Af' (n + ,)/2n f 


p x + a y 


= A f( n_1 )/ 2n 
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f 3 = 4f(l - f) 
f =2(1-20 

(n+l)/2n 

f (* + ^)x =(l/2)"Af+ AAf 

(n + 1 )/2n 

f (or?+P*)£ = - (1/2)R A f +* Af 

or., y' + p = («• ? I- B ? ) f/(( l/2)rr f + P f ) 

1+y ' p = (n* + P°)f/((l/2hf + Pff 
(3n + l)/2n 

Ay"=MR P T f 4}/2)^r + 8ff 


Direct substitution shows that indeed 

(3n-l)/2 

(y*') n =(n-l) Vy' + pXl+y 1 *) 

For the special case n = 2, the parameter t is related to curvature k by the relation 
k = cos t. Arc length s is related to the parameter t by the relation s = (cos t) 1 , 
and strain energy E(t) at any point is related to curvature by the relation E = K. 

Each of the previous cases is a special case of a more general approach. In 
general, we can let 



34 


We observe then that 


a x - B y =[An/(n-l)] F 
px + nr, y =[An/(n-l)] G 
(*"+ P^x = [An/(n-l)] F +fl G) 
(a^+pp); =[ An /(h-l)] (-B F + n G) 

ory' + p =( a '+r)G/( (1 'F+pG) 

(3n-l)/(l-n) 

•••••• , _ • /9 ' " ' ' 

GF-FG =[n/(n-l)]G F 

1 + y"= («•* + p")(F% G")/(a F+ pdf 

„„ , <3n-l)/(l-n) . , 

y“ = -("’+f)GF /A("F + 0 G) 

F%df =<rF ,n/(, - n> 


and direct substitution verifies that (y") n = (n-1) (ay'+ P)(l + y' 2 ) 


(3n-1)/2 


This permits an infinite variety of parametric representations. We need only 
choose F arbitrarily and let G be given by 


BOUNDARY CONDITIONS 


The problem of matching the solutions of the foregoing differential equations 
was to have been the major emphasis of the second phase of this project. Con- 
sequently, as of the writing of this summary report, very little attention has been 
devoted to the boundary condition problem. 

Basically, the problem is as follows: in each of the intervals (x x^) 
n = 1,2, , N we have a solution of the form f (x, y, n - n , p n , y^, fi n ) = 0. 

Consequently, we have 4N constants of integration to be determined. With these 
constants, we have the conditions that the solutions pass through the specified 
points (2N conditions), that y 1 be continuous at each joint (N-l conditions) and 
that y" be continuous at each joint (N-l conditions). In addition, we impose 
additional condition at each end (at X q and at x^). Consequently, we have 4N 
specified conditions to determine the 4N constants of integration. While, in theory, 
this is sufficient, the practical situation is that this is a formidable problem indeed. 
For example, in several of the situations the solutions involve elliptic functions 
whose invariants are dependent upon the constants of integration. In these cases, 
functional values cannot be determined except as transcendental functions of the 
constants of integration. So, in general, the problem reduces to one of solving 4N 
non-linear equations in 4N unknowns. Detailed study of individual cases might 
well result in effecient individual procedures for attacking this problem. But, as of 
this report, no such detailed study has been attempted. 
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It is observed that, in each case, the differential equation provides an 
additional route for investigation of half the constants of integration (the and 8^). 

In this process, however, one must temporarily introduce additional unknown 
quantities (the unknown values of the slope and the second derivative at each joint). 

The following very general approach may indicate a possible process for determining 
the constants of integration. Of the eight differential equations we have considered, 
the six which result from approximate integrands are of the form 

y"=(ct y' + p ) f(y') x€ (x , x ) , n = 1,2, , N 

n 1 n n-i n 

If we assume that y" has no zeroes, except possibly at the joints, and introduce the 
quantities 

m n = y'(x ) n =0, 1, 2, , N 

n 

and 

M n = y" (x ) n = 0, 1 , 2, , N 

n 

we obtain the "parametric solutions" (for n =1,2, , N) 

y' 

x-x = S l |"( a n U + )f(*j) du 

• n-j n J 

m'* L J 

n-i 

and 
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M n-1 - ^ a n M n + n = l,2, — , N 

M; - (a n M i P n )f(M ) n = 1,2, N 

n n n n n 

and two end conditions. This, in theory, would permit one to determine the 4N+2 
unknown quantities if the values of s n were supplied in advance. Some obvious 
simplifications are possible but it is not presently possible to state whether or not 
this general process would be either effecient or fruitful . We can only conclude 
that the boundary condition problem must be given considerable attention before 
it is possible to make a valid conclusion as to the utility of our solutions for 
applications oriented users. 
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APPENDIX A: SPECIAL CASES 


In our discussion of the eight differential equations given on page 5, we have 
generally excluded the cases where nr. f = 0. The purpose for this exclusion was that, 
generally, four arbitrary constants were required in each interval in order to satisfy 
the imposed boundary conditions. In the case that a particular constant is zero by 
viture of boundary conditions, then the solution will appropriately reduce to fit 
the boundary conditions. In many physical applications, continuity of y" at the 
joints is essential (e.g. in automotive or ship design problems, the mathematical 
curves are reproduced by mechanical means. The milling (or drafting) machine 
will not operate in such a way as to allow the second derivative (and hence the 
currvature) to vary in a discontinuous manner. Such variance would require 
machinery to achieve instantaneous acceleration.) However, in many curve fitting 
applications, continuity of second derivatives may not be an essential feature. This 
fact has led us to briefly examine the results of setting iQf= 0 in the differential 
equations in question. 
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The Special Cases a = 0 

The cases which arise by putting a = 0 in the differential equations have a 

special interest since they can also be interpreted as the results of minimization 

of the strain energy integral when continuity of y" is ignored. We saw earlier 

' b 

that the Euler-Lagrange equations for extremizing the integral $ T(y , /y ,, )dx / when 

a 

y" cPyii = 2co, was f P = ay' + p. If we put y's u in the integral, we seek to extremize 

f b 

a '3 ^(i^u'Jdx. The Euler-Lagrange equation is then cp= u' fpyi - p . But u' ip, u' 3 2 <p # So, 

we have rp - /3 as the differencial equation. This is the same result that is obtained by 
putting a = 0 in original case. 

In this case, the equations from the six integrands given on page 3 become, 
respectively, 

1. y"Vp n 

2. y"~~ p (1 +5y |S /2) 

3. y"° = P n /( 1 -5y'*/2) 

4. y"°=P (l+5y ,S /4f 

n 

5. y" S = S /(l-5y ,S /4f 

n 

6. y"° = 0 n (i + y ,3 f /3 

we note that, in each case, P n > 0. For simplicity of notation, we delete the sub- 
script n and change the notation of the arbitrary constant by, respectively, taking 
f n = Pn = 2P = /S, P n =5PV 2 , p n = <4p/5 f, P n = (5 p/4 f , and f n -p\ This 
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gives the six equations 


1 . 

2 . 

3. 

4. 

5. 

6 . 


y" = f 


y“'‘= 8 n (a%y* ? ); a 2 =2/5 
y" = f /(a - y* ) ; a' = 2/5 
y" = P (a" + y*‘ ) ; a' = 4/5 
y ,,P= P7(a p -y ,? f ; a 7 =4/5 
y ,,? = pWV* 


As in previous studies, we assume that y" has no zeroes in the interval in question. 

Consequently, we may assume that 6 and u" have the same sign and f 4 0. 

2 2 . 

The equation y" = p. 

With our assumptions, this gives y" = f* and y is a second degree polynomial in x. 
Note that there are three arbitrary constants to be determined 
The equation y" = p (a y 1 ") 

In this case we have y" = f (a'-i- y' ) "whose solution 
is p y + * = a cosh ( x + y ) (where a > 0) . 

The equation y"' = £ /(a -y'") 

In this case we have y"(a -y 1 ' ) = f* . Now, since y" =dy'/dx and y” = y'dy'/dy 

we have 

9 x i- ^ $ (a - u ) ^ du 


and 


Py + Y = (a-u^udu 
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where ft and V are constants of integration and the parameters is, in fact, y'. 

Now, in each integral we change variables by letting u = a sin 9 (where 
_ ir/2< 6 < 5r/2). We obtain then 

"J3 x + ft = a 3 $cos 9 d9 

and 


py+ v 



sin 9 d9 


We can readily integrate these equations to obtain 


|8x+ ft = (a /2) (9 + sin 9 cos 9 ) 


Px + Y = -(a /3) cos 9 

It is easily verified that this parametric solution does indeed satisfy the differential 
equation. It is also possible to eliminate the parameter to obtain a functional relation 
between y and x. Since a > 0 and cos 9 > we must have ( 8 y + Y) ^ 0. Then 

a cos 9 = -[3(f?y +Y)] l/s , 

9 = s cos 1 lj-3/a" )( p y + Y)J ^ , 
and a sin 9 = s j a*- [3(j3 y + y)] P/ ^ | 

where s = * 1 =sgn(y'). Consequently, 

2s( P x + M=° COS ’ [ (-3/cT)( Py + v)] 1 ^ 

- [3( e y + V) ] ^ (p y + } l/s 
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The Equation y"' = ft ? (a + y r ) 


In this case, we have y" = fl (a °+ y'~). Again, since y" = dy'/dx 
and y" = y' dy'/dx, we have 

P x + fi = J(a + u ) du = (l/a) tan *(u/a) 

and 


p y 1 Y * = 5 ( a? u ) udu ~ 0/2) ^ n ( a + u ) 
where and V* are constants of integration and the parameteru is, in fact, y*. We 
note that these equations can be written 


and 


Thus, 


u/a = tan |a( fix + f>j| 
(1/2) l n (1 + u~/a ) = p x + y 
f y + T = (1/2) £n [sec a( fx + K)| 


or, 

P y t y = - 4 n cos a( p x + ^) 

The Equation y" = ft /(a - y 1 ) 

In this case, we have y"(a°- y'*) = p. Again, since y" = dy'/dx and y" = y’ 
dy'/dx, we have 


<• a •-» q 

p x + f> - J (a' - u )du = au - u /? 
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and 


fy+Y = (a - u") du = (a“u / 2 ) - u/4 

where 6 and Y are constants of integration and the parameter u is, in fact, y'. 

Now, from the second equation, we find that tT = a i [ a “4 (Py + Y)] */ 2 

l/rr 

However # u 2 -a S «r* 0 so that u ? = a - [a 4 - 4(^y + Y)j . Now, we find by combining the 

original two equations that a ( fl x -i- y )/2 - ( f y + V) -- u /l 2. But, 
u* = 2a 4 - 4( Py + y) - 2a 2 [a -4(Py i y)J . So, we have 

u^ = a~ - [a 4 - 4( p y + y)] ^ 

and 

3u (px+S) = a 4 + 4(p y + ^-a'' [a 4 - 4( p y + y)J 1 / s 

If the second equation is squared, and then u is replaced by the first expression, 
we obtain a functional relation between x and y. 

The Equation y" - [3 (1 + y 1 ") ^ 

In this case we have y" = (3(1 I y 1 ') . By inspection, we see that equation 

is satisfied when ever 

y r = [(-P y/2) +v] "*-i 

We may assume, with no loss of generality that 0 <r (-py/2) + Y 1 . Now, we put 

-l A 

(- (3y/2) + y = cn (s,k) where k = 2 “ . With this substitution we obtain 

(p dx/ds)" = 2 cn 4 (s) 

so that 
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fix + * = + 2 V ^J crT(s) ds 


So, 


+ * = + 2^ Je(s) - s/2] 


where E(s) is the fundamental elliptic integral of the second kind. 

There does not appear to be a method of eliminating the parameter s so 
as to achieve a functional relation between x and y. 
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BOUNDARY CONDITIONS 


In the preceeding special cases, the boundary condition problem is quite 
different from that of the general case. In these cases, there are three constants 
per interval. Generally, two of the constants (say f, and v) can be eliminated 
by requiring the solution to fit through the prescribed joints. This leaves, for 
N + 1 intervals, N + 1 f's to be determined. The requirement of continuous 
first derivatives at the joints will provide N conditions on these N + 1 constants 
and there remains only one additional condition that can be imposed. In this 
case then, we cannot provide the usual two end conditions. Since the approx- 
imated integral can be integrated and expressed in terms of the R's one possible 
approach is to take a fixed value for the (approximated) total strain energy as 
the final condition. This might be of particular interest in computer graphics 
where one could generate, and compare, an infinite family of curves by varying 
this one variable. 

As mentioned earlier, there is every indication that the boundary condition 
problem will be a formidable problem. As of the writing of this report, very little 
attention has been devoted to this aspect of the general problem. 
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